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Abstract.  A  procedure  developed  for  free  molecule  modeling  of  plume  backflow  from  a  STAR™  30BP  solid  rocket 
motor  is  presented  for  work  perfonned  in  support  of  the  Comet  Nucleus  Tour  spacecraft  mishap  investigation.  Good 
general  agreement  is  established  with  DSMC  flowfield  results,  with  interesting  deviations  developing  as  the  plume 
backflow  approaches  the  spacecraft  surfaces  closely,  providing  insights  regarding  characteristics  of  the  surface  Knudsen 
layer.  Also,  investigation  of  related  free  expansion  results  indicate  significant  discrepancies  exist  between  the  rarefied 
techniques  and  the  continuum  results  from  which  their  starting  surfaces  were  created.  The  nature  of  these  differences 
suggests  that  convective  fluxes  to  CONTOUR  may  have  been  much  higher  than  the  rarefied  analyses  indicated. 


INTRODUCTION 

On  15  August  2002,  the  NASA/Johns  Hopkins  University  Applied  Physics  Laboratory  (APL)  Comet  Nucleus 
Tour  (CONTOUR)  spacecraft  was  destroyed  during  conclusion  of  an  ATK  STAR™  30BP  solid  rocket  motor 
(SRM)  orbit  transfer  burn1.  A  NASA  Mishap  Investigation  Board  (MIB)  was  established  one  week  after  loss  of 
contact,  and  under  its  direction,  a  number  of  failure  modes  were  investigated  and  revisited  in  greater  detail1.  Without 
telemetry  coverage  or  direct  optical  observation  to  guide  them,  the  MIB  was  unable  to  determine  the  cause  of  failure 
indubitably1.  However,  they  eventually  established  the  most  probable  cause  through  a  remarkable  sequence  of 
events,  beginning  with  overheating  of  CONTOUR  by  backflow  from  the  SRM  exhaust  plume1.  The  SRM  nozzle 
was  nearly  immersed  within  the  spacecraft  (Fig.  1),  allowing  nearby  surfaces,  some  of  which  were  bare  AE,  to  be 
exposed  to  intolerably  high  heat  fluxes.  Such  “....overheating  may  have  caused  substantial  material  weakening  and 
structural  degradation,  which  could  have  led  to  catastrophic  dynamic  instability”1. 


FIGURE  1.  Representation  of  CONTOUR  spacecraft,  including  SRM  and  SPAF  [2], 

The  MIB  decision  was  guided  by  sets  of  analyses  describing  convective  and  radiative  surface  heating  due  to 
SRM  operations2'5.  In  the  final  round,  one  organization  used  CFD  to  model  SRM  product  expansion  from  the 
nozzle  to  near  field  plume  development  in  order  to  provide  a  starting  surface  for  rarefied  techniques  and  to  establish 
convective  backflow  heating  estimates  from  the  continuum  standpoint3.  The  starting  surface  provided  input 
conditions  for  DSMC/DS2V  computations  as  well  as  a  free  molecule  technique  for  the  same  problem4'5. 
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This  paper  describes  the  free  molecule  (FM)  analysis5,  an  approach  which  has  previously  fared  well  in 
comparisons  with  a  variety  of  DSMC  simulations  and  experiments6,7.  Good  general  agreement  is  established  with 
the  two-dimensional  DSMC/DS2V  flowfield,  with  interesting  deviations  developing  as  spacecraft  surfaces  are 
approached.  These  differences  illustrate  the  Knudsen  layer’s  influence  in  suppressing  the  excursion  of  surface- 
reflected  molecules  away  from  the  body.  Finally,  comparisons  of  related  free  expansions  indicate  significant 
discrepancies  exist  in  the  rarefied  techniques  and  the  continuum  results  from  which  their  starting  surfaces  were 
created.  The  nature  of  these  differences  suggests  that  convective  fluxes  to  CONTOUR  may  have  been  much  higher 
than  the  rarefied  analyses  indicated. 


MODEL  FORMULATION 


A  transient  solution  of  the  collisionless  Boltzmann  equation  was  developed6'7  to  describe  the  molecular 
distribution  fix,t)  for  flow  from  a  point  source  step  function  Qu  where 
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dt  dx 
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Q\  represents  directed  flow  from  a  Lambertian,  thermal  velocity  distribution  superimposed  on  convective  exit 
velocity  ue  for  mass  rate  m  ,  / 6  =  l/ *J2RT  ,  and  speed  ratio  s  =  P  ue  .  Normal  n  represents  the  orientation  of  a  local 
starting  surface  element,  and  v  •  n  emphasizes  the  imposed  directional  constraint.  Generally,  ne  is  not  aligned  with 
h ,  with  the  angle  between  the  two  defined  by  </>e.  Local  angle  (/  is  measured  between  variable  position  x  (distance  r, 
experiencing  local  velocity  v)  and  n ,  and  angle  6  is  measured  between  «e  and  x. 

The  particular  solution  of  Eq.  (1)  is  found  using  approaches  outlined  by  Bird*  and  Narasimha9.  The  density  field 
generated  for  a  step  function  in  m  is  given  by6'7 


p{x,t)  = 


/3  m  cost/)  w2  2 

- 

Al  nr~ 


(2) 


where  z  =  a  -  w,  a  =  [3  r/t ,  and  w  =  .v  cos  8  .  Solving  for  successive  velocity  moments,  one  finds  mass  flux  ®  , 
incident  normal  momentum  flux  p  \ ,  and  incident  translational  energy  flux  gTR  : 
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Eq.  4  becomes  normal  momentum  flux  p  when  multiplied  by  the  dot  product  r  makes  with  the  local  normal  of 
the  impinged  surface.  Equations  (2)-(4)  may  be  combined  to  obtain  expressions  for  velocity  v,  translational 
temperature  FTR,  and  internal  energy  flux  c/|Nt  for  polyatomic  molecules  having  specific  heat  ratio  y. 
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Reflected  quantities  such  as  normal  momentum  flux  and  energy  are  found  by  setting  s  =  0,  letting «  represent  the 
local  surface  element,  and  assuming  the  mass  flux  to  a  surface  element  is  conserved8.  The  expansion  around  the 
nozzle  lip  into  the  backflow  region  may  be  portrayed  by  using  the  near  field  solution  of  a  continuum  code  to  create  a 
starting  surface  for  a  network  of  FM  point  sources  having  locally  varying  conditions. 


PROCEDURE 

CONTOUR  consisted  of  an  octagonal  structure  2.3  m  across  its  apothems  on  the  dustshield  end  (Fig.  1).  The 
SRM  was  located  within  a  load-bearing  cylinder  featuring  the  spacecraft  payload  attachment  fitting  (SPAF).  The 
SPAF  cylinder  extended  13  cm  beyond  the  dustshield,  and  the  SRM  nozzle  extended  just  another  6  cm  beyond  an 
annulus  of  protective  insulating  blankets.  Unprotected  metal  blanket  buttons  and  an  aluminum  low-gain  antenna 
(LGA)  were  exposed  to  the  plume  heating.  Surface  temperature  was  assumed  constant  at  333  K. 

In  early  MIB  analyses,  protruding  assemblies  were  neglected,  and  the  overall  geometry  was  modeled  as  three 
stacked,  concentric  cylinders.  Two  cases  were  modeled  in  the  FM  investigations;  one  including  the  body,  and  a  free 
expansion.  In  the  former,  logic  for  one  set  of  diffuse  surface  scattering  was  added  to  observe  its  perturbation  on  the 
free  expansion,  in  addition  to  the  perturbation  due  to  physical  blocking  effects  from  the  structure.  The  STAR™ 
30BP  motor  consumes  9.2  kg/s  of  material  at  a  specific  impulse  of  292  s  for  about  50  s  10,  and  SRM  gaseous  fuel 
products  have  been  listed  in  Table  1  below. 


Table  1.  STAR™  30BP  SRM  Model  Exhaust  Composition  at  2450  K10 


Species 

h2o 

h2 

CO 

O 

O 

HCf 

n2 

ce 

H+ 

Gaseous  Mass  Fraction 

0.156 

0.028 

0.300 

0.052 

0.329 

0.130 

0.005 

<0.001 

Specific  Heat  Ratio  y 

1.18 

1.30 

1.29 

1.16 

1.30 

1.29 

1.65 

1.67 

Actual  SRM  exhaust  would  also  contain  liquid  and  solid  aluminum  compounds.  These  should  be  confined  to  a 
relatively  narrow  cone  once  ejected,  and  were  ignored  for  the  convective  backflow  analyses. 

Exhaust  composition  is  important  because  as  the  plume  expands,  a  species  separation  effect  occurs,  allowing 
lightweight  constituents  to  spread  more  diffusely  while  heavier  ones  remain  more  tightly  focused  along  the  plume 
axis6'7'12.  Such  effects  become  increasingly  important  as  one  moves  to  higher  angles  off  the  plume  axis. 


RESULTS  &  DISCUSSION 


FIGURE  2.  Comparison  of  steady  DSMC/DS2V  (a)  and  current  FM  results  (b),  logarithmic  number  density  contour  maps. 
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Figure  2  presents  a  comparison  of  steady  number  density  («)  contour  maps  for  a  radial  cross-section  on  the  same 
logarithmic  scale.  There  is  a  nice  level  of  agreement  overall,  considering  that  for  the  DSMC  case  the  nozzle  lip 
region  was  recalculated  beginning  slightly  ahead  of  the  exit  plane4,  and  more  detail  has  been  included  in  the  nozzle 
and  SPAF  regions  there. 

On  the  other  hand,  not  all  differences  between  the  two  results  can  be  explained  by  lower  input  resolution  for  the 
FM  case.  As  spacecraft  surfaces  are  approached,  it  is  clear  that  n  increases  in  the  DSMC  case,  while  it  begins  to 
drop  for  the  FM  case.  Although  the  author  did  not  have  access  to  raw  DSMC  simulation  output  data,  careful 
measurements  were  made  from  Fig.  2a  above  at  radial  distances  of  0.4,  0.6,  0.9,  and  1.2  m  from  the  SRM  central 
axis  and  compared  to  the  FM  results  (Fig.  3). 


FIGURE  3.  Comparison  of  steady  flowfield  number  density  variations  for  DSMC/DS2V  and  current  FM  results  at 
(a)  0.4  m,  (b)  0.6  m,  (c)  0.9  m,  and  (d)  1.2  m.  CONTOUR  spacecraft  surfaces  located  near  the  vertical  axis  of  each  plot 
(left-hand  side),  starting  surface  located  along  right-hand  side.  FM  free  expansion  results  given  by  lower  curve. 

Although  it  was  sometimes  too  difficult  to  resolve  the  precise  axial  location  of  individual  DSMC  contour  levels 
near  the  body  (Figs.  3b  &  3c),  they  qualitatively  follow  the  behavior  depicted  in  Fig.  3a.  The  rise  in  density  for 
DSMC  results  is  associated  with  the  radial  development  of  the  Knudsen  layer  across  the  SPAF  closeout  ring  and  the 
dustshield.  Without  the  interaction  between  surface -reflected  molecules  and  those  approaching  directly  from  the 
plume  starting  surface,  the  FM  solution  shows  no  difference  between  this  case  and  the  free  expansion  (no  body)  at 
the  body  locations. 

This  absence  of  FM  surface-reflected  molecules  right  at  the  body  location  is  a  consequence  of  using  a  diffuse, 
cosine  distribution  for  the  scattered  ensemble,  which  assumes  all  such  molecules  must  have  a  component  of  velocity 
moving  away  from  the  body.  Local  density  is  then  distributed  farther  away  from  the  body,  giving  rise  to  higher 
density  levels  relative  to  DSMC.  These  balance  out  the  deficits  observed  contiguous  to  the  surfaces.  Researchers 
generally  refer  to  the  portion  of  the  DSMC  profile  in  Fig.  3a  where  density  increases  near  the  body  as  the  edge  of  the 
Knudsen  layer.  Without  a  Knudsen  layer,  FM  results  show  an  enhancement  of  density  further  afield. 

The  development  of  a  Knudsen  layer  in  the  DSMC  solution  gives  rise  to  relative  variations  in  heat  flux 
distributions  versus  the  FM  results  (Fig.  4  below).  Both  solutions  seem  to  rise  from  near  zero  to  0.28  W/cm2 
midway  across  the  SPAF  closeout  ring,  but  then  the  influence  of  the  lip  in  the  DSMC  case  causes  fluxes  to  increase 
further,  reaching  about  0.45  W/cm2  before  the  lip  is  reached.  This  geometric  detail  was  neglected  for  FM  runs. 

Moving  radially  along  the  dustshield,  both  solutions  reach  a  peak  net  heat  transfer  of  about  0.1  W/cm2  before 
declining  to  about  0.05  W/cm2  at  its  outer  edge.  However,  dustshield  peak  values  occur  at  a  lower  radius  for  the  FM 
case.  It  may  be  possible  to  explain  this  shift  either  by  suggesting  higher  DSMC  cell  resolution  may  be  necessary,  or 
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that  some  sort  of  local  radial  convective  process  is  delaying  heat  flux  inputs  to  the  dustshield  through  the  Knudsen 
layer. 


FIGURE  4.  Comparison  of  surface  convective  heat  transfer  balance  versus  running  length  away  from  SRM  nozzle  lip. 
Net  heat  transfer  presented  in  yellow  for  DSMC  case  (a,  middle),  purple  for  FM  case  (b,  second  from  top). 


One  alarming  discovery  made  during  this  investigation  regarded  free  expansion  solutions  created  by  the 
MOC/boundary-layer/CFD  approach  and  the  current  FM  technique  (Fig.  5).  Somehow,  the  anticipated  exponential 
decay  of  mass  flux  with  angle14  was  not  observed  for  FM  results  using  the  CFD  starting  surface,  even  though  the 
CFD  free  expansion  from  which  the  starting  surface  was  created  follows  this  relationship  itself. 
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FIGURE  5.  Comparison  of  STAR™  30BP  mass  flux  angular  dependency  for  steady  free  expansions, 
CFD  results  and  FM  results  using  CFD  starting  surface,  performed  ~  1  m  off  centerline. 


No  free  expansion  calculations  were  performed  by  the  DSMC  investigator,  but  the  good  agreement  between  FM 
and  DSMC  density  contour  maps  presented  in  Figs.  2-4  suggests  a  DSMC  free  expansion  would  behave  in  a  fashion 
similar  to  the  FM  result,  since  fluxes  of  order  higher  than  density  typically  behave  roughly  like  a  product  of  density 
with  some  suitable  power  of  ultimate  velocity.  Discrepancies  of  10-20x  in  mass  flux  would  certainly  be  evident  in 
density  comparisons  such  as  Figs.  2  &  3. 

Based  on  these  observations,  it  seems  reasonable  to  assume  that  the  depressed  mass  flux  versus  angle 
relationship  at  high  angles  is  due  to  a  deficit  in  density  in  precisely  the  same  region  where  SRM  backflow  impinges 
on  CONTOUR.  Such  a  deficit  would  yield  proportional  effects  on  incident  heat  fluxes,  which  is  in  keeping  with 
early  FM  work  for  the  CONTOUR  MIB  using  an  analytically-based  starting  surface  (Fig.  5)5.  Based  upon  these 
observations,  one  must  consider  that  actual  heat  fluxes  impinging  on  CONTOUR  in  the  backflow  regime  may  have 
been  10-20x  higher  than  what  had  otherwise  been  predicted. 

Even  using  the  heating  rates  presented  in  Fig.  4,  other  follow-on  analyses2  indicated  that  structural  integrity 
could  have  weakened  enough  during  the  SRM  burn  to  cause  CONTOUR  to  fail  by  precipitating  dynamic  instability. 
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CONCLUDING  REMARKS 


A  free  molecule  rocket  plume  expansion  model  was  used  to  estimate  incident  heat  fluxes  in  the  backflow  regime 
of  an  SRM  on  an  actual  spacecraft.  FM  flowfield  density  results  generally  compared  favorably  with  DSMC 
calculations,  except  as  the  body  was  approached  closely.  In  the  absence  of  a  Knudsen  layer,  free  molecule  results 
suggest  that  flowfield  behavior  near  such  surfaces  is  dominated  by  the  gas-surface  interaction  model  employed,  and 
that  the  flowfield  is  thus  affected  over  a  greater  region  than  what  is  generally  identified  with  the  Knudsen  layer. 

Attempts  to  verify  an  expected  mass  flux  versus  angle  relationship  revealed  significant  discrepancies  exist 
between  FM  and  DSMC  solutions  versus  the  continuum  results  from  which  their  starting  surfaces  were  created. 
Further  investigation  did  not  determine  the  cause  behind  this  density  deficit,  but  it  is  considered  possible  that  actual 
heat  fluxes  impinging  on  CONTOUR’S  SPAF  closeout  ring  and  dustshield  may  have  been  even  higher  than  the 
significant  levels  presented. 
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